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The red clump is found to be split into two components along several sightlines 
toward the Galactic Bulge. This split is detected with high significance toward 
the areas (—3.5 < / < 1,6 < —5) and {l,b) = (0, +5.2), i.e., along the Bulge 
minor axis and at least 5 degrees off the plane. The fainter (hereafter "main") 
component is the one that more closely follows the distance-longitude relation of 
the Bulge red clump. The main component is ~0.5 magnitudes fainter than the 
secondary component and with an overall approximately equal population. For 
sightlines further from the plane, the difference in brightness increases, and more 
stars are found in the secondary component than in the main component. The 
two components have very nearly equal {V — I) color. 
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1. Introduction 

Investigations of the red clump (RC) have ha d a strong imp a ct on studie s of the Galactic 
Bulge. In their analysis of OGLE-I photometry, IStanek et al.l (119941 . 119971 ) found an east- 
west asymmetry in the dereddened apparent magnitudes of the RC, suggesting a bar oriented 
at 20-30° relative to the line o f sight toward the Galac tic center (GC), with the near side 
in the first Galactic quadrant. iRattenburv et al.l (120071) obtaine d similar results using the 
much larger OGL E-II survey, as did Babusiaux fc Gilmord (120051) using near -IR photometry 
along the plane. iBenjamin et al.l ( 20051 ) and ICabrera-Lavers et al.l ( 20071 ) independently 
found evidence for a long-bar component with a length of ~4.5 kpc and a viewing angle 
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of ~45°. ICabrera-Lavers et al.l (120081 ) argued further that a triaxial spheroid component 
coexists with the bar, and dominates in the range |/| 5, 10°. Further ev i dence of multiple 
dynamical components are the works of lAlardI (120011 ) and lNishiyama et al.l (|2005[ ). who found 
evidence that the stellar population toward the range (|/| ^ 4, 6 < 1) is distributed differently 
than that of the main bar. 



The OGLE-III survey has the potential to yield some of the most powerful constraints 
yet on the major questions: Are there one or two Bulge populations? Where does each 
dominate? What are their viewing angles and characteristic lengths? With ~96 deg^ of 
observations in V and /, the OGLE-III survey is the largest photometric database of the 
Galactic Bulge that consistently penetrates deep enough to isolate the RC in the coordinate 
range (|/| < 10, 2 < |6| < 7). In the course of fitting for the position of the RC we found 
that there frequently would appear to be two RCs in a sightline, with the best-fit magnitude 
being near the mean of that of the two components. We concluded that this effect must be 
better understood before full modelling of the Bulge can be undertaken. 

In this Letter, we describe the double red clump of the Bulge in the context of our 
observations: where it is detected, and the relative properties of the two populations in 
magnitude, color, dispersion and number counts. 



2. Data 



OGLE-III observations were obtained from the 1.3 meter Warsaw Telescope, located at 
the Las Campanas Observatory in Chile. The camera has eight 2048x4096 detectors, with 
a combined field of view of 0.6° x 0.6° yielding a scale of approximately 0.26"/pixel. We 
use observations from a mosaic of 267 fields directed toward the Galactic Bulge, which are 
almost entirely within the range —10 < / < 10 and 2 < |6| < 7. The photometric coverage is 
shown in Figured] Of the 267 fields, 37 are at northern latitudes. More detailed descriptions 
of the instruni entat ion, photomet r ic red uctions and astrometric calibrations are available in 
Udalskil J2OO3I ) and ludalski et al.N2008h . 



Method 



We fit for the position of the RC in a multi-step process. We first estimate the position 
of the RC by eye, by looking at sample locations within every OGLE-III field. This is used 
as a guess for the following step. At each of ~92,000 windows spaced approximately 2' apart, 
a geometric centroid in (/, V — I) is calculated based on stars that are within 0.5 magnitudes 



CD 
Zi 



o 



o 



5 - 



- 



5 - 



-3- 



"1 — I — r 



"1 — I — I — r 



"1 — r 



<ff 






"1 1 



# 






Wf 



A.. 



"1 — I — I — r 



1 — r 



CT„p < 0.270 
0.270 < CTrc < 0.289 




0.306 
0.326 
351 
393 



< 
< 
< 
< 



RC 
RC 
RC 
RC 



< 0. 

< " 

< 0, 



326 
393 






."t^ 










J L 



10 



5 0-5 

Galactic Longitude 



-10 



Fig. 1. — The best-fit width of the RC obtained assuming a single Gaussian. The distribution 
shown as a color-map of septiles of equal areas. The contours of the areas used in this study 
are dehneated by the black rectangles - any OGLE-field with a sub-field centered within 
those contours is used. Areas where it was difficult to fit for the RC due to factors such as 
differential reddening are left out as white space. 



of the guessed clump location on the color magnitude diagram (CMD). Each location is 
inspected by eye to verify the fit, and the process is repeated until every inspection window 
is found to have a visually acceptable RC centroid. It was in this step, with the large number 
of direct visual inspections required, that we first noticed that a number of sightlines had 
multiple RCs in their CMDs. 

A more detailed analysis is then required, to delineate the peak of the RC population 
from the peak of the combined RC-fRG (Red Giant) population, and to estimate properties 
such as the number density of RC stars and their dispersion in magnitude within a given 
sightline. To account for these factors, we select stars for fitting that satisfy the following 
color- magnitude cut: 



0.2 <{V-I)-{V~ I)rc,g < 0.2 
-1.5 < / - IncG < 1.0 



(1) 
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where the subscript (RC,G) refers to the guessed RC location inferred from the previous 
process of geometric centroiding. The RC+RG population is then fit with the equation: 



N{m) = Aexp B{m- rriRc) H — 1= exp —^ , (2) 

2ai 



'RC 



Nrc 

"Itxorc 

where Nrc is the number of RC stars within a window, vtirc is their estimated brightness 
peak, grc is their dispersion in brightness, and A and B fit an exponential profile to the 
magnitude distribution of RG stars, corresponding to a power-law fit to the flux. The best-fit 
is obtained by minimizing x^ over 50 histogram bins in the magnitude range of Equation 

2 V^ {Nmodel — Nobs) /„n 

There is a degeneracy that manifests itself as the width of the RC aRc gets larger. The 
distribution of stars can then be well-fit by a range of parameter values involving larger 
aRc, larger Nrc, and lower B, i.e., it becomes harder to distinguish between the RC and 
the RG branch as the RC gets wider. As such, we impose a Gaussian prior, (/x, cr)/- = 
(0.3652, 0.0269), on the fraction /j of stars in a given sightline that are RC stars divided by 
the sum of Nrc and the RG stars that satisfy the criteria of Equation ([T]). 

This procedure proves adequate for large swaths of the Galactic Bulge, examples are 
shown in Figure [2l However it is not adequate for regions where the RC is split into two. 



3.1. The Split Red Clump 

Our first computational evidence of the existence of a split RC is the distribution of aRc 
as a function of longitude and latitude, plotted in Figure [TJ There is an east- west asymmetry 
in the southern Galactic hemisphere. In the coordinate range (1 < / < 5,6 < —5), aRc is 
typically ~0.3 magnitudes, whereas for the coordinate range (—3 < / < 1,6 < —5), aRc 
exceeds 0.4 magnitudes. On closer inspection, we observed is that there are two "bumps" 
toward many of these sightlines. The Gaussian fit is then spurious, with its "peak" falling in 
between the two "real" peaks, and its width broadened to partially encompass both bumps. 
Examples are shown in Figure [3l This effect may also be present in other parts of the sky, 
but if so the two clumps are likely separated by less than 0.4 magnitudes, rendering it more 
difficult to distinguish them. 

We decided to focus on those OGLE-HI fields where the double RC would be most 
prominent: BLG125, 126, 127, 128, 134, 135, 136, 137, 143, 144, 145, 146, 151, 152, 153, 
154, 159, 160, 161, 162, 167, 168, 169, 170, 176, 177, 178, 186, 355, 357; approximately 
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corresponding to the areas (—3.5 < / < 1,6 < —5.0) and (—2 < / < 2,5 < 6). A full list of 
OGLE-III Bulge fields and their coordinates can be found on the OGLE webpagqj. 



In each window, we fit the RC+RG population to the following equation: 



N{m) = Aexp 



B{m-mRc,i] 



N, 
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Within the color-magnitude selection box: 



exp 



'^'^ RCA 



2-, 



N, 



RC,2 



2'Ka 



exp 



RC,2 



0.2 <{V-I)~{V- I)nc,G 
-1.5 < I -Ircg < 1-5 



[m - mRc,2) 

•^^ RC,2 

(4) 



(5) 
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The angular radius used is the smallest one such that ~1200 stars meet the conditions of 
Equation ([5]), and the coordinates selected for inspection are placed 3' apart within every 
OGLE-III subfield - they will be closer if they are at the edges of neighboring fields. Degen- 
eracies and correlations remained a concern in early trial runs, and thus we imposed three 
priors on the fit. The fraction of stars fi that are RC stars (of either clump) and meet the 
conditions of Equation ([5]) was given the Gaussian prior of (/i,cr)/. = (0.25,0.05); Gaussian 
priors were also used to discourage solutions with a clump dispersion a^c < 0.1 mag, or a 
brightness peak either at least 1 magnitude brighter or 0.5 magnitudes fainter than the peak 
obtained from a single Gaussian fit. 



3.2. Comparing V — I color 

We compare the V — I color of the two RC by inspecting 621 high- confidence sightlines 
for which A^^ > 8 when a second Gaussian is added to the fit, the difference in magnitude 
is 0.4 < A/ < 0.6, and the smaller RC has at least one quarter the population of the larger 
RC. The effect of these selection criteria is shown in Figure HI In each sightline, a box is set 
with an area of (0.3 x 0.3) mag^ in color-magnitude space, centered on each of the two RCs 
in the window. The difference in color, {V — I)bright — {V — I) faints had a mean of 0.0034 
magnitudes and a dispersion of 0.012 magnitudes. It is thus consistent with zero. Values 
ranged between —0.028 mag and +0.036 mag, which is roughly the 3o" error range, where a 
is the error in the mean determined determined from the typical dispersion in l^ — J color 
of the stars in a box (~0.08 mag) and the typical number of stars in a box (~120). 



^http://ogle.astrouw.edu.pl/^ogle/ 
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3.3. Properties of the Two Red Clumps as a Function of Reddening, 

Longitude and Latitude 

A least-squares fit is computed for the magnitude as a function of reddening, longitude 
and absolute latitude on the same 621 selected sightlines, using the following formalism: 

/ = /o + i?,*i?(y-/) + ^*Z + ^*(|6|-5), (6) 

wher e Rt = At/E(V — I) is ratio of total to selective extinction, E{V — I) assumes {V — I)o = 
1.08 ( JBensby et al.l 120101 ) . and -£ and ^ parametrize any change in distance modulus to 



the source populations due to their spatial distribution. The values obtained are: 

Imain = 14.735 + 1.168 * E{V -I)- 0.0136 * / + 0.018 * (|6| - 5) 
Isecondary = 14.259 + 1.168 * E{V -I)- 0.0136 * / - 0.010 * (|6| - 5) (7) 

We fixed Rj and ^ to be the same for both populations as they only differed at the ~2.0cr 
and ~1.0cr level respectively. By contrast, the latitude terms, differ by ~6.0cr. The ~2cr 
detection of a difference in Rj may be due to systematic effects such as the presence of dust 
between the two populations. 

4. Results 

We observe the following: 

• The main cl ump better traces the dereddened Bulge RC magnitude known from previ- 



ous studies (IStanek et al.lll997l : iRattenbury et al.l 120071 ). though it is ~0.1 magnitudes 
fainter. 

The secondary RC is about 0.5 magnitudes brighter than the main RC. The difference 
in brightness increases further from the plane. 

The two RC are approximately equally populated, with the secondary RC becoming 
more and more populated relative to the main RC further from the plane. 

The two RCs have approximately equal {V — I) color. 

The main RC becomes narrower for sightlines further from the plane, whereas the 
secondary RC becomes broader. 
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• The same trends are observed in the northern hemisphere, but they are less statistically 
significant since there are 15x fewer data. 

• The secondary RC becomes more and more populated relative to the main RC as 
longitude increases. 

• The difference in brightness Imam — Secondary, and the ratio of the widths cr_Rc,mam '■ 
o'RC, secondary, show uo Statistically significant trends with longitude. 

• Within the OGLE-III Bulge database, the effect is visible in the coordinate range 

(—3.5 ^ ^ ^ 1, 5 < |6|). If it is present elsewhere, it is less pronounced. 

In Figure [5|, we plot the distribution of the ratio of the populations of the two RCs, 
NRCmain/NRC secondary, and the difference in brightness Imain - ^secondary as functions of 
southern Galactic Latitude . 



Discussion and Conclusion 



It is not clear what this double clump feature could be. iGallart fc Bertellil (Il998l ) sum- 
marized many of the pit-stops on the red giant branch, such as the RGB-bump and the 
AGB-bump. The AGB-bump for example, would seem to fit the bill as it is suggested to be 
of a similar color to the red clump and about 1 magnitude brighter. However, why would it 
be visible in these specific parts of the sky? It is also not expected to be equally populated 
to the RC. 

An investigation of the RC across the entire OGLE-III Galactic Bulge data set is in 
progress (Nataf et al. 2010, in preparation). It may be found that the expected brightness 
at the locations used for this study is in between those of the two RCs. If that is the case, 
one possible explanation is that one RC belongs to the Galactic bar and the other to the 
spheroid component of the Bulge, and that due to geometrical effects their separation is 
more obvious for (—4 < / < 1, |6| > 5). 

It is not even certain that this double RC effect is localized - it is only certain that its 
degree of detectability is localized. We have deviations between the two clumps of between 
0.4 and 0.6 magnitudes. If elsewhere on the sky they were separated by 0.4 magnitudes or 
less, they might be indistinguishable from a single broad clump. Philosophically, the signal 
is not from two peaks but from a trough between the two peaks. This trough is suppressed 
by statistical noise and by the underlying RG branch, which is a monotonically increasing 
function of magnitude. The intermediary distribution between two peaks with a trough and 



a single peak, i.e., that of a flat distribution over a range of ~0.4 magnitudes, could also be 
produced by a boxy Bulge. As such, if it were seen it could not be argued to be necessarily 
due to the presence of the two RCs moving closer together. 



The color constraint is a powerful one. iGirardi fc Salarid (120011) investigat ed the behav- 
ior of the RC as a function of age and metallicity using the iGirardi et al.l ( l2000l ) evolutionary 
tracks and isochrones. They found that for metallicities ranging from Z = 0.0004 through 
to Z = 0.03, the (y — I) color levels off after about 3 Gyr. The brighter RC could, for 
example, be both of a higher metallicity than the regular RC and younger than 2-3 Gyr - 
the two effects could then cancel out and the difference in /-band would be consistent with 
the observed value of 0.5 magnitudes. If this secondary RC is older than ~3 Gyr, the two 
populations would need approximately equal metallicity, and the secondary RC would need 
to be closer to explain the magnitude difference. For example, for respective ages of the 
populations of 3 and 10 Gyr, the secondary RC would need to be about 1.0 kpc closer. For 
ages of 7 and 10 Gyr, it would ne ed to be about 1.5 kpc closer. We note that the model 
results in IGirardi fc SalarisI (120011 ) assume a scaled-solar metallicity - a broader range of 
choices may be available if [a/Fe] is allowed to differ between the two RCs. 

A captured and cannibalized dwarf galaxy is another possibility. Since its stars would 
occupy a specific set of orbits, this could explain why this feature is seen in some parts of 
the sky and not o thers. The current nearest galaxy to the MW is the Sagittarius stream 
( llbata et al.lll995l ): it is at an earlier stage of digestion than this hypothetical stream would 
be. Radial velocities (RV) would likely be significantly different between the two populations 
in such a scenario and as such RV measurements would provide an effective constraint on the 
nature of this population. Table [T]lists 5,085 RC stars including their OGLE ID, astrometric 
coordinate, photometry and whether or not they are of the main or secondary component. 
Proper motions (PM) are an alternative option to the same ends. PM measu rements are 



not a s precise as those done with RV. The proper motion catalog for OGLE-II ( ISumi et al. 



20041 ) includes two fields that overlap with the sightlines used in this study, BUL_SC6 and 
BUL_SC7. A proper reduction catalog of OGLE-III (Poleski et al. 2011, in preparation) will 
yield measurements for tens of thousands of RC stars in the reg ions of interest . Currently, we 
find only one mention of the Bulge double RC in the literature: IZoccalil (120101 ) briefly reports 
on a double RC finding along the bulge minor axis and at distances in excess of ~700 pc 
removed from the pl ane. An early analysis of Spitzer-GLIMPSE mid-IR photometry closer 
to and in the plane (JBenjamin et al.l 120031 ) also shows the presence of some multiple peaks 
(Benjamin et al. 2010, in preparation). 

We bring up two other "double RC" findings in the l iterature for comparison. The first 
being the globular cluster Terzan 5 (JFerraro et al.l |2009| ) and the other in the Sagittarius 
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streain (jCorrenti et al.ll2010l ). Neither case is directly analogous to our own. The two RCs 
in Terzan 5 are separated not just in magnitude, but in color as well, with the fainter RC 
being about ~0.15 magnitudes bluer in (J — K). This lends itself to an explanation of 
multiple stellar populations with differ ent ages and iron cont ent, or with similar ages and a 
different iron as well as helium content (JD'Antona et al.ll2010l ). The Sagittarius dwarf galaxy 
on the other hand has a more natural ge ometric explanation t o its double RCs - sightlines 
intersecting with multiple wraps of stars. ICorrenti et al.l ( I2OIOI ) detect multiple double RCs 
often separated by 0.5 magnitudes or greater. 



Moving forward, one way to ascertain if this double clump is widespread is to obtain 
multi-color photo r netry as well as spectra in regions where the two clumps are distinct. 
Salaris fc Girardil ( 120021 ). analyzing the results of their simulations in i^-band, found that 
the brightness of the clump increases with increasing metallicity, whereas V and /-band 
brightness decrease with increasing metallicity. If distinctions such as an enhanced metallicity 
appear, they could be used to confirm or rule out the possibility that this double RC is more 
common in the Bulge than just in the regions of interest where we verified its presence. The 
implications would be wid espread, as the width of the RC, aRc is used to constrain models 



of th e Bulge morphology ( IStanek et al.l 119971 : iRattenbury et al.l l2007l : ICabrera-Lavers et al. 
20081 ). If the width is due to the presence of multiple stellar populations, rather than 



a dispersion due to a range of distances, then the thickness of the Bulge obtained from 
RC studies would be seriously overestimated and the viewing angle of the bar would be 
underestimated. Of equal importance, variations in apparent magnitude would then often 
be due to the varying relative size of the two populations, rather than the a change in 
distance modulus due to the orientation of the bar. 

DMN and AG were partially supported by the NSF grant AST-0757888. The OGLE 
project has received funding from the European Research Council under the European 
Community's Seventh Framework Programme (FP7/2007-2013) / ERG grant agreement no. 
246678 to AU. 
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further study. A briefer version of the table is shown below. Selected stars had a {V — I) 
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Fig. 2. — Color-Magnitude diagrams and brightness distribution with best-fit single clump 
model for two locations. Left: (/, 6) = (1.0,-3.9), with paramaters {rnRc,(^RC,NRc,) = 
(15.177,0.273,288.6). Fitting a second Gaussian improves x^ from 54.3 to 50.7, i.e. 3 
additional d.o.f. yields A^^ = 3.6. Right: (/,6) = (—4.0,-3.3), with clump paramaters 
{jriRc 1 (^RC 1 Nrc i) = (16.085,0.360,310.1). A^^ = 3.1 for a second Gaussian. 
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Fig. 3. — Color-Magnitude diagrams and brightness distribution with best-fit double 
clump model for two locations. Left: (/, 6) = (0.27,-5.77), the Ax^ = 26.7 when a sec- 
ond Gaussian is fit, the parameters are: (^_rc,1' '^-RC.i? ^RC,i;''^J?c,25 o"_rc,2, ^^0,2) = 
(14.90,0.18,166.42,15.42,0.10,98.51,). Right: (/, 6) = (-0.28,5.76), the 

Ax^ = 14.4. The parameters are: 
(15.42, 0.19, 147.46, 15.91, 0.10, 94.23, ). 
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Fig. 4. — TOP: The dereddened main and secondary RC are respectively plotted in black 
and red as a function of longitude for every siglitline. Since a double-RC is not prominent 
in every siglitline a few sequences are seen, including a possible red giant branch bump at 
I ~ 15.2. BOTTOM: The plotted points are those for which a second Gaussian yields a 
Ax^ > 8 and the smaller RC with at least one quarter the population of the larger RC. 
Those for which the difference in magnitude is 0.4 < A/ < 0.6 are plotted as filled circles. 
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Fig. 5. — The magnitude difference between the bright and faint clump, and the ratio of the 
two populations as a function of latitude. We plotted points for those 621 high- confidence 
sightlines described previously. The medians are shown as red points. 



